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Abstract 

Until recently, the heat capacity of GaN had only been measured for polycrystalline powder 
samples. Semiempirical as well as first-principles calculations have appeared within the past few 
years. We present in this article measurements of the heat capacity of hexagonal single crystals 
of GaN in the 20-1400K temperature range. We find that our data deviate significantly from the 
literature values for polycrystalline materials. The dependence of the heat capacity on the isotopic 
mass has also been investigated recently for monatomic crystals such as diamond, silicon, and 
germanium. Multi-atomic crystals are expected to exhibit a different dependence of these heat 
capacities on the masses of each of the isotopes present. These effects have not been investigated 
in the past. We also present first-principles calculations of the dependence of the heat capacities of 
GaN, as a canonical binary material, on each of the Ga and N masses. We show that they are indeed 
different, as expected from the fact that the Ga mass affects mainly the acoustic, that of N the 
optic phonons. It is hoped that these calculations will encourage experimental measurements of the 
dependence of the heat capacity on isotopic masses in binary and more complex semiconductors. 

PACS numbers: 65.40. +g; 63.20.Dj 
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I. INTRODUCTION 



The investigation of the heat capacity of crystals is an old topic of condensed matter 
physics with which illustrious names were early associate d 1 ' 2 ' 3 . Knowledge of the heat ca- 
pacity of a substance not only provides essential insight into its vibrational properties but 
is also mandatory for many applications. 

Two famous limiting cases are correctly predicted by the standard elastic continuum 
theory 3 -: At high temperatures, the constant- volume heat capacity, C v , tends to the Petit 
and Dulong limit, e. g. ~ 49.9 J/mol K in binary materials such as GaN^. At sufficiently 
low temperatures, C v is proportional to T 3 (Ref. (^)). At intermediate temperatures, 
however, the temperature dependence of C v is governed by the details of vibrations of the 
atoms and for a long time could only be determined from experiments. In recent years, the 
temperature-dependence of the specific heat of many systems has been calculated in the 
harmonic approximation. The calculations ranged from semiempirical, based on inelastic 
neutron scattering or X-ray data, to first-principles, based on total energy calculations. 

For example, for GaN, a very topical semiconductor these days, the temperature depen- 
dence of the specific heat capacity has been calculated, in the harmonic approximation, on 
the basis of either experimental or based on first-principles phonon frequencies. The phonon 
frequencies were obtained semiempirically, from inelastic neutron 5 - or x-ray scattering data£ 
or from first-principles, based on electronic calculations of total energies^. 

In the past 15 years a variety of semiconductor crystals with different isotopic composi- 
tions have been grown 8 . Using these crystals, the effects of the isotopic mass on the phonon 
spectra^, in particular the related anharmonic effects, have been investigated. Of consid- 
erable interest has also been the dependence of the electronic spectra on isotopic mass, 
which is effected by the electron-phonon interaction^. The influence of the isotopic masses 
on thermo-mechanical properties, including the thermal conductivit y 11 ' 12 ' 13 ' 14 the thermal 
expansion 15 and the specific hea t 1 ^' 1 ?' 18 have also been investigated. Monatomic materials 
such as diamond, Si, and Ge have only one average isotopic mass that can be varied (effects 
of fluctuations in the isotopic mass will not be discussed here). GaP and GaAs are binary 
compounds but their anions have only one stable isotope, hence only the effects of varying 
the cation isotope mass can be investigated (even so, very interesting effects have been seen 
in the Raman phonons of GaP (see Ref. Il9j)). 
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GaN is an ideal material to investigate isotope effects: Ga possesses two stable isotopes, 
( 69 Ga and 71 Ga), whereas N has also two such isotopes, 14 N and 15 N. An investigation of the 
electronic energy gap of GaN and its dependence on the isotopic mass of N was published 
recently^. Corresponding measurements in which the isotopic mass of gallium is varied 
were, at that time, not possible because of the lack of suitable isotopically modified samples: 
they are being performed now. 

Isotopically modified samples of GaN large enough to measure the dependence of the 
heat capacity on the masses of Ga and N are still not available. In this paper we present 
first-principles results concerning this problem, in the hope that they will encourage crystal 
growth with isotopically pure Ga and N. Our theoretical results are displayed, in the Debye 
manner, as C v /T 3 , for the four possible isotope combinations, compared with the calculations 
for GaN with the natural isotopic abundances. The maximum in C v /T 3 found at about 40 
K increases when any of the isotopic masses is increased, this effect for gallium being larger 
than for nitrogen. The effect of the latter becomes larger than that of the former at higher 
temperatures. We also present the calculated logarithmic derivative of C v /T 3 vs. each of the 
two masses: the effect of the mass of nitrogen is shown to be displaced to higher temperatures 
with respect to that of the mass of gallium. This can be qualitatively attributed to the fact 
that the gallium mass affects mainly the acoustic phonons of GaN whereas that of N affects 
mainly the optical phonons. The temperature dependence of the logarithmic derivatives of 
C v /T 3 is compared with the results obtained for monatomic crystals, in which only one mass 
can be varied. 

Due to the difficulties involved in growing large enough single crystals, the experimental 
results available thus far for GaN have been obtained for polycrystalline powder samples 
onl y 21 i 22 i 23 i 24 . In this article, we present new measurements of the heat capacity in the range 
20-1400 K of single-crystals of GaN grown by the high-pressure, high-temperature technique. 
Our new heat capacity data are in good agreement with the results of the first-principles 
calculations. However, particularly above room temperature, the heat capacities previously 
obtained on powder samples, and commonly quoted in the literature as standard values, 
deviate significantly from our data and first-principles theory. We therefore suggest that 
the standard values for the heat capacity of GaN should be reconsidered and replaced. 
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II. EXPERIMENTAL 



GaN single crystals were grown by a self-seeding process at pressures of nitrogen of about 
1.5 GPa and temperatures of 1800 K from a nitrogen solution in a droplet of galliu m 25 i 26 i 27 . 
The crystals used in this experiment were of up to 5 mm 2 lateral size and thickness of 40- 
300 fi. They have a good crystallographic quality (x-ray diffraction rocking curves of 30-100 
arcsec). They were highly conductive (n > 5xl0 19 cm -3 ) due to unintentionally introduced 
oxygen donors but showed a very low concentration of dislocation which did not exceed 10 2 
cm -2 . After removing the crystals from the high pressure growth chamber they were etched 
in aqua regia which can result in a traces of Ga20"3 on the etched surfaces of crystal plates. 

The heat capacities in the temperature range 15 K < T < 400 K were measured using 
a commercial Physical Property Measurements System calorimeter (Quantum Design, 6325 
Lusk Boulevard, San Diego, CA.) employing the relaxation method. To thermally anchor 
the crystals to the sample holder platform, a minute amount of Apiezon N grease was 
used. The total heat capacity of the platform and grease was determined in a separate 
run and subtracted from the measurements for the GaN crystals. The heat capacities in 
the temperature range 300 K < T < 750 K were measured with a Perkin-Elmer Pyris 1 
Differential Scanning Calorimeter (Perkin Elemer, 45 William Street, Wellesley, MA 02481- 
4078 U.S.A.) in a Pt crucible. A total of three independent runs were merged and fitted with 
a polynomial. The deviation of the experimental points from this polynomial was always 
less than ±1.5%. The maximum observed difference between data points for the different 
runs was ~2%. In the temperature range 350 K < T < 1400 K we used a Netzsch 404C 
Pegasus DSC calorimeter (NETZSCH-Geratebau GmbH, Wittelsbacher Str. 42, D-95100 
Selb, Germany) with the sample contained in a Pt crucible equipped with alumina liners. 
For the DSC measurements several crystals were combined to obtain samples with masses 
between 20 and 60 mg. During the measurement the DSC calorimeters were flushed with 
N2 gas. The calibration of the DSC calorimeters was done with a sapphire crystal. 

III. FIRST-PRINCIPLES CALCULATIONS 

The theoretical results described below were obtained from self-consistent, first-principles 
calculations based on local density-functional theory with the SIESTA code 28,29 , as de- 
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scribed in Ref. |30j The exchange-correlation potential of Ceperley- Alder— , parameterized 
by Perdew and Zunger— , was used. Norm- conserving pseudopotentials in the Kleinman- 
Bylander forrn 3 ^ were used to remove the core regions from the calculations. The basis sets 
for the valence states were linear combinations of numerical atomic orbitals (double-zeta 
with a set of d orbitals for Ga). The host crystal was represented by a (wurtzite) Ga36N 36 
periodic supercell. The calculated lattice constants and bulk modulus were found to be a = 
3.200 A, c/a = 1.629, and B = 211 GPa, respectively, which compare well with the literature 
values 34 (3.189 A, 1.607, and 202 GPa, respectively). 

The matrix elements of the (harmonic) dynamical matrices were extracted at T = K 
from the derivatives of the density matrix relative to nuclear coordinates using the perturba- 
tive approach developed by Baroni et a/.— and Gonze et alM and implemented into SIESTA 
by Pruneda et al. 37 . The dynamical matrix was then evaluated at 120 q points along high- 
symmetry directions of the Brillouin zone. This procedure generated some 2.6xl0 4 normal 
mode frequencies. The phonon density of states (pDoS) shown in Fig. Qwas obtained by 
smoothing the normal mode frequencies with Gaussians with a width of 1 cm -1 . 

In the harmonic approximation, the Helmholtz free energy is given by 3 ^ 

,{T) = k B T \n{smh(hu/2k B T)}g{u;)du; (1) 
Jo 

where k B is the Boltzmann constant. Here, the pDoS g{uj) is normalized so that J g{uj)duj = 
3N, where N is the number of atoms. Note that Fmb.iT = 0) is the total zero-point energy. 
Once F V ib. is calculated, the vibrational entropy and specific heat at constant volume are 
given by 



(^) C v = - T (^\ . (2) 



dT J v ~ u ' V ~9T 2 
In Fig. EJ we display C v /T 3 calculated in the 20-60 K range for six isotopic compositions, 
including natural Ga ( natural nitrogen can be assumed to be pure 14 N for our purposes). 
This figure clearly illustrates that the effect of changing the mass of Ga is considerably larger 
than that of changing the mass of N in the temperature region under consideration. This 
is because the vibrations of nitrogen take place at much higher frequency than those of Ga 
(see FigllJ. Figure 01 shows the logarithmic derivatives of C v /T 3 vs. the masses of either 
the Ga or the N atoms, as well as vs. temperature calculated from the data of Fig. |21 This 
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plot was particularly useful for monatomic crystals because it could be related to either the 
measurecU&i£ or calculated temperature dependence of specific heat by the equation^: 

dln(C v /T 3 ) = d\n(C v /T 3 ) 3 

dlnM d\nT 2" U 

We have also plotted in Fig. El a curve obtained with Eq. El using the data of Fig. HI In 
comparing this plot with the calculated K values of both derivatives, we realize that the 
results of Eq. Q correspond roughly to the sum of the two separate derivatives. This is 
understood by considering that in the monatomic case the masses of both atoms in the unit 
cell are varied when taking the mass derivative, whereas in the binary case we vary only one 
at a time. In the limit (T — > ) Eq. (jHJ) yields a derivative equal to 3/2. The logarithmic 
derivatives vs. mass in the binary case (e.g. GaN) are then given by: 



dlnC v 3 M Ga 



dlnM Ga 2M Ga + M N 
din a, 3 Mm 



1.25, (4) 



0.25. (5) 



dlnM N 2M Ga + M N 

The sum of both derivatives in Eqs. ()4|5)1 equals 3/2, as predicted for a monatomic 
crystal. This trend can already be seen in the derivatives shown in Fig. El for 20 K. In 
this figure, the gallium derivative remains larger than the nitrogen derivative up to about 
160 K. At this temperature the contribution of the gallium to C v has already flattened out 
considerably, while tending to the Petit and Dulong limit, whereas that of nitrogen has 
begun to rise up and takes over. 

Notice that below 80 K the derivative vs. Mn increases down to 20 K, the lowest tem- 
perature we could reliably reach in our calculations. This can be attributed to a dominant 
contribution of long- wavelength acoustic phonons which include vibrations of Ga as well as 
of N (the appropriate effective mass is M Ga + Mn). The maximum in the Mn derivative 
seen in Fig. El at ~ 220 K corresponds to that found for the M Ga derivative at 30 K. It 
signals the rising contribution of the nitrogen related acoustic phonons. 
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IV. EXPERIMENTAL RESULTS AND DISCUSSION 



Figure E| displays our low temperature molar heat capacity data of crystalline GaN as 
a Debye plot (C p /T 3 vs. T), together with the results of the first-principles calculations 
in the temperature range 20 K < T < 300K. For these temperatures, thermal expansion 
contributions to the heat capacity are negligible, and to a very good approximation, C p ~ C v 
(cf. the discussion below). 

The agreement of our results with the first-principles calculations is very good, small 
deviations (< 2.5%) are visible below the maximum at ~ 40 K. The magnitude of the Debye 
maximum is well reproduced in the experimental data, T max is shifted by 1.5 K to higher 
temperatures from the value predicted by the calculations. For comparison we also show the 
data by Koshchenko et ali^ obtained from measurements of a polycrystalline sample. Below 
~ 200 K Koshchenko's data, taken so far as the standard values for the low temperature 
heat capacity of GaN, deviate significantly from the first-principles results and the present 
experimental data. They completely fail to reproduce the maximum at ~ 40 K, which is 
due to deviations from a T 3 power law behavior at low temperatures. We presume that 
impurities may be the reason for this failure. Above 150 K, Koshchenko's results agree 
fairly well with our data and the first-principles results. The heat capacities determined by 
Leitner et al. for temperatures 200 K < T < 1300 K (not shown in Fig. 0} consistently lie 
above our data. Up to 300 K the difference does not exceed 2%. 

We discuss next our heat capacity results in the the temperature range 300 K < T < 
1400 K where the difference between const ant- volume and constant-pressure specific heats 
becomes increasingly important. The constant- volume and constant-pressure specific heats 
are related by^ 

C P (T) = C V (T) + E(T)-T, (6) 

with 

E(T)=a 2 v (T)-B-V moh (7) 

where a v (T) is the temperature dependent volume coefficient of thermal expansion, B the 
bulk modulus and Vmoi the molar volume at T = K. The linear thermal expansion coeffi- 
cients, a a {T) and a c (T), have been measured by x-ray diffraction over a wide temperature 
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range (100 K < T < 1500 K) by Wang and Reeber— . The volume thermal expansion co- 
efficient, a v (T) = 2 a a (T)+a c (T), amounts to 12.29xl0~ 6 K _1 at 300 K and it increases 
to 15.57xl0~ 6 K _1 at 1500 K. Above 200 K the overall temperature dependence of the vol- 
ume thermal expansion coefficient can be well described by a v (T) = 15.59(4) xlO -6 K -1 -(l 
- 0.25(1) xl0 5 K 2 /T 2 + 0.50(8) xl0 9 K 4 /T 4 ). The overall agreement of this fit with the ex- 
perimental results is better than ±1%. 

Figure El shows our measured molar heat capacities, C p , of crystalline GaN in comparison 
with the data available in the literature for polycrystalline GaN. The data by Leitner et al. 22 
and Yamaguchi et alM, obtained on powder samples and commonly quoted as standard 
values for the specific heat of GaN at high temperatures, and in particular the data by Chen 
et al—^ deviate considerably from our results. Above 800 K, our data taken on a pristine 
sample, reveal a broad hump centered at about 900 K. This anomaly is no longer seen in a 
subsequent second run carried out on the same sample (see the discussion below). Above 
1100 K, the results of both runs agree again within experimental error. 

The agreement of our high temperature heat capacity data with C p = C V + E(T) ■ T, with 
C v obtained from first principles by Sanati and Estreicher^, is within the expectations for 
such calculations. To calculate the thermal expansion contribution E(T) according to Eq. 



©, the bulk modulus B = 202.4 GPa (Ref. |34j), the molar volume V moi = 13.76xl0~ 6 m 



,3 



and the polynomial fit of the temperature dependent a v (T) described in detail above were 
used. 

The solid (black) line in Figure El was obtained by least-squares fitting all our data (for 
T > 300 K, excluding the anomaly around 900 K seen in the first run) simultaneously to a 
power series expansion of the Debye integral, including terms oc 1/T 2 and oc 1/T 4 as well 
as a term C • T to account for the thermal expansion. For simplicity we have assumed a 



temperature independent Debye temperature 9^ (Ref. |42j). 



1 6 2 1 9 4 

CJT) = 2 x 3R (1 ^ + rs)+C7-T, (8) 

p\ ) v 20 T 2 560 T 4 ; ' K ' 

with R being the gas constant. 

Our measured molar heat capacity C v in the temperature range 300 K < T < 1400 K 
can be best described with 0^ = 863(3) K. This Debye temperature correlates very well 
with the 0£>(T) obtained from the measurements below room temperature and is consistent 
with 0£>(T) calculated by Nipko et a/.— from inelastic neutron scattering data, (see inset in 



Fig. El). 

The fit with Eq.fjHJ gives C = 1.03±0.01 mJ/mol K 2 , somewhat larger than the value of 
~ 0.7 mJ/mol K 2 expected using the thermal expansion coefficient for T — > oo given above. 

There is an anomaly centered at ~ 900 K, but starting already at ~ 700 K, in the heat 
capacity data collected in the first run up to 1400 K. This feature was no longer seen in a 
subsequent run carried out on the identical sample. Additionally, we detected a small weight 
loss starting at ~ 600 K which finally amounts to ~ 0.05% at 1400 K. The heat capacity 
anomaly and the weight loss indicate sublimation, most likely of the Ga-suboxide species, 
Ga20. Sublimation of Ga20 in high-vacuum has been found to commence at 500 °C^i. 
Since the DSC measurements were carried out in N 2 atmosphere Ga 2 could result from a 
reduction of adherent traces of Ga203 detected after the freshly grown crystals have been 
etched with aqua regia (see above). 

In summary, we determined the heat capacity of GaN single crystals in the tempera- 
ture range 20 K < T < 1400 K and found significant deviations from the values given in 
the literature for polycrystalline material. Good agreement of our experimental data with 
first principles calculations is observed over the full temperature range. We make detailed 
predictions as to how the heat capacity will vary if the isotope masses of Ga and N are 
changed. We hope that these interesting results will soon be confirmed by measurements on 
isotopically pure GaN crystals. 
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FIG. 1: Calculated phonon density of states (pDoS) of GaN taken from Ref. [fl. The band below 
400 K corresponds basically to Ga vibrations (except in the acoustic region below ~ 100 K where 
both atoms vibrate). The optical band, above 500 K, corresponds to nitrogen vibrations. 
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FIG. 3: (color online) Calculated logarithmic derivatives dln(C 1) /T 3 )/(ilnT + 3/2 (right-hand- side 
of Eq. © (solid black line) compared to the left-hand-side d ln(C v /T 3 )/d In M with M = M Ga 
(short-dash line, blue) and M = Mn (long-dash line, red). Note that the sum of the two dashed 
lines is roughly equal to the solid black line. 
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FIG. 4: (color online) Low temperature experimental molar heat capacity C p /T 3 of wurtzite GaN 
(circles). The (red) solid line represents the results of the first-principles calculations for C v /T s 
(C v ~ C p ). The (blue) dashed line represents the data obtained by Koshchenko et al. on a 
poly crystalline sampled. The inset shows the region of the Debye maximum in an enlarged scale. 
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FIG. 5: (color online) The experimental heat capacities C p of wurtzite GaN crystals compared with 
theoretical results. The solid (black) line represents the best fit of Eq. (JSJ) to our data combining 
all our data (except the anomaly around 600 K, see text). The solid (green) line represents the 
first-principles C v data by Sanati et al. (Ref. |7j). The solid (red) line shows C p as calculated 
according to Eqs. Q6I7|) using the literature values for the bulk modulus and the molar volume 
and the expansion coefficients for the volume thermal expansion, a v (T), as given in the text. The 
high temperature data collected with the Netzsch calorimeter (•, o) show an anomaly at ~ 900 K, 
most likely due to the sublimation of traces of Ga20 (see text). A concomitant small weight loss 
has been observed in a TGA trace in the respective temperature region. The dotted (black), the 
dashed (black and the dash-dotted lines show literature data by Leitner et al. (Ref. 22), Itagaki 



et al. (Ref. 



23f ) and Chen et al. (Ref. |2J), respectively. The inset displays the Debye temperatures 



corresponding to our data according to the Debye law (20 K < T < 400 K), in comparison with 
Nipko's results**-. To fit the high temperature data with Eq. (jSJ) a temperature-independent Debye 
temperature 0oo, indicated by the horizontal bar, has been assumed. 
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